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Abstract 

The model of weak measurements is apphed to various problems, related to the 
time problem in quantum mechanics. The review and generalization of the theoretical 
analysis of the time problem in quantum mechanics based on the concept of weak 
measurements are presented. A question of the time interval the system spends in 
the specified state, when the final state of the system is given, is raised. Using 
the concept of weak measurements the expression for such time is obtained. The 
results are applied to the tunneling problem. A procedure for the calculation of 
the asymptotic tunneling and reflection times is proposed. Examples for (5-form 
and rectangular barrier illustrate the obtained results. Using the concept of weak 
measurements the arrival time probability distribution is defined by analogy with the 
classical mechanics. The proposed procedure is suitable to the free particles and to 
particles subjected to an external potential, as well. It is shown that such an approach 
imposes an inherent limitation to the accuracy of the arrival time definition. 



1 Introduction 



The time plays a special role in quantum mechanics. Unlike other observables, time remains 
a classical variable. It cannot be simply quantized because, as it is well known, the self- 
adjoint operator of time does not exist for the bounded Hamiltonians. The problems related 
to time also arise from the fact that in quantum mechanics many quantities cannot have 
definite values simultaneously. The absence of the time operator makes this problem even 
more complicated. However, in practice the time is often important for an experimenter. 
If quantum mechanics can correctly describe the outcomes of the experiments, it must also 
give the method for the calculation of the time the particle spends in some region. 

The most-known problem of time in quantum mechanics is the so-called "tunneling 
time problem". Tunneling phenomena are inherent in numerous quantum systems, ranging 
from an atom and condensed matter to quantum fields. There have been many attempts 
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to define a physical time for tunneling processes, since this question has been raised by 
MacCoU [1] in 1932. This question is still the subject of much controversy, since numerous 
theories contradict each other in their predictions for "the tunneling time" . Some of these 
theories predict that the tunneling process is faster than light, whereas the others state 
that it should be subluminal. This subject has been covered in a number of reviews (Hauge 
and St0vneng [2], 1989; Olkholovsky and Recami [3], 1992; Landauer and Martin [4], 1994 
and Chiao and Steinberg [5], 1997). The fact that there is a time related to the tunneling 
process has been observed experimentally [6-15]. However, the results of the experiments 
are ambiguous. 

Many problems with time in quantum mechanics arise from the noncommutativity of 
the operators. The noncommutativity of the operators in quantum mechanics can be cir- 
cumvented by using the concept of weak measurements. The concept of weak measurement 
was proposed by Aharonov, Albert and Vaidman [16-21]. Such an approach has several 
advantages. It gives, in principle, the procedure for measuring the physical quantity. Sec- 
ond, since in the classical mechanics all quantities can have definite values simultaneously, 
weak measurements give the correct classical limit. The concept of weak measurements 
has been already applied to the time problem in quantum mechanics [22-24]. 

The time in classical mechanics describes not a single state of the system but the process 
of the evolution. This property is an essential concept of the time. We speak about the time 
belonging to a certain evolution of the system. If the measurement of the time disturbs 
the evolution we cannot attribute this measured duration to the undisturbed evolution. 
Therefore, we should require that the measurement of the time does not disturb the motion 
of the system. This means that the interaction of the system with the measuring device 
must be asymptotically weak. In quantum mechanics this means that we cannot use the 
strong measurements described by the von-Neumann's projection postulate. We have to 
use the weak measurements of Aharonov, Albert and Vaidman [16-21], instead. 

We proceed as follows: In Sec.|21we present the model of the weak measurements. Sec. El 
presents the time on condition that the system is in the given final state. In Sec. |3J our 
formalism is applied to the tunneling time problem. In Sec. El the weak measurement of 
the quantum arrival time distribution is presented. Section (HI summarizes our findings. 

2 The concept of weak measurements 

In this section we present the concept of weak measurement, proposed by Aharonov, Albert 
and Vaidman [16-21]. We measure quantity represented by the operator A. 

We have the detector in the initial state |$). For a weak measurement to provide the 
meaningful information the measurements must be performed on an ensemble of identical 
systems. It is supposed that each system with its own detector is prepared in the same 
initial state. After measurement the readings of the detectors are collected and averaged. 

Our model consists of the system S under consideration and of the detector D. The 
total Hamiltonian is 

H = Hs + Ht, + Hi (1) 

where Hs and are the Hamiltonians of the system and detector, respectively. We 
take the operator describing the interaction between the particle and the detector of the 
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form [22, 23, 25-28] 

Hi = XqA, (2) 

where A characterizes the strength of the interaction between the system and detector. The 
small parameter A ensures the undisturbance of the system's evolution. The measurement 
duration is r. In this section we assume that the interaction strength A and the time r are 
small. The operator q acts in the Hilbert space of the detector. We require the spectrum 
of the operator q to be continuous. For simplicity, we can consider this operator to be the 
coordinate of the detector. The momentum which conjugate to q is Pq. 

The interaction operator ((21) only slightly differs from the one used by Aharonov, Albert 
and Vaidman [17]. The similar interaction operator has been considered by von Neumann 
[29] and has been widely used in the strong measurement models (e.g., [28,30-34] and 
many others). 

Hamiltonian Q represents a constant force acting on the detector. This force results 
in the change of momentum of the detector. From the classical point of view, the change 
of the momentum is proportional to the force acting on the detector. Since interaction 
strength A and the duration of the measurement r are small, the average {A) should not 
change significantly during the measurement. The action of the Hamiltonian Q results in 
the small change of the mean detector momentum {pq) — {pq)o = —Xt{A), where {pq)o = 
($(0)|pq|$(0)) is the mean momentum of the detector at the beginning of the measurement 
and (pq) = (<l'(r)|pg|$(r)) is the mean momentum of the detector after the measurement. 
Therefore, in analogy to Ref. [17], we define the "weak value" of the average (A), 

(A) = MiLziM. (3) 

AT 

At the moment t = the density matrix of the whole system is /3(0) = ps(0) Pd(0), 
where ps(0) is the density matrix of the system and Pd(0) = |$)($| is the density ma- 
trix of the detector. After the interaction the density matrix of the detector is p-oit) = 

Trs (ps(0) (S> I $)($!) f/^(t)| where U{t) is the evolution operator. Later, for simphc- 

ity we shall neglect the Hamiltonian of the detector. Then, the evolution operator in the 
first-order approximation is [22] 

f/(t)^f/s(t)(l + ^£i(ti)dti^ (4) 

where Us(t) is the evolution operator of the unperturbed system and A{t) = Ul(t)AUs{t). 
From Eq. Q we obtain that the weak value coincides with the usual average {A) = 
Tr{ips(0)}. 

The influence of the weak measurement on the evolution of the measured system can be 
made arbitrary small using the small parameter A. Therefore, after the interaction of the 
measured system with the detector we can try to measure the second observable B using, 
as usual, the strong measurement. As far as our model gives the correct result for the value 
of A averaged over the entire ensemble of the systems, now we can try to take the average 
only over the subensemble of the systems with the given value of the quantity B. We 
measure the momenta Pq of each measuring device after the interaction with the system. 
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Subsequently, we perform the final, postselection measurement of B on the systems of our 
ensemble. Then we collect the outcomes Pg only of the systems which have a given value 
oiB. 

The joint probability that the system has the given value of B and the detector has 
the momentum pq(t) at the time moment t is W{B,pq]t) = Tt {\B){B\pq){pq\p{t)}, where 
|pq) is the eigenfunction of the momentum operator pq. In quantum mechanics the proba- 
bility that two quantities simultaneously have definite values does not always exist. If the 
joint probability does not exist then the concept of the conditional probability is meaning- 
less. However, in our case operators B and |j5q)(j9q| act in different spaces and commute, 
therefore, the probability W{B,pq;t) exists. 

Let us define the conditional probability, i.e., the probability that the momentum of 
the detector is p^ provided that the system has the given value of B. This probability is 
given according to Bayes theorem as 

where W{B]t) = TT{\B){B\p{t)} is the probability that the system has the given value 
of B. The average momentum of the detector on condition that the system has the given 
value of B is 

(PqW) = j P^W{p^;t\B)dp^. (6) 

IFtoiu. Eqs. (jHI) and (jH)), in the first-order approximation we obtain the mean value of 
A on condition that the system has the given value of B (see for analogy Ref. [22]) 

{Ah - ^^^{\BmA+A\Bm 



^ihimfW) ^^^^^^^^ " ^''^^^^^^ ( [i^^^^i'^] ) • 



ih{B\ps\B) 

If the commutator \B){B\,A in Eq. ((Tj) is not zero then, even in the limit of the very 

weak measurement, the measured value depends on the particular detector. This fact 
means that in such a case we cannot obtain the definite value. Moreover, the coefficient 
{{q){Pq) — Re{qpq)) may be zero for the specific initial state of the detector, e.g., for the 
Gaussian distribution of the coordinate q and momentum pq. 

3 The definition of time under condition that the sys- 
tem is in the given final state 

The most-known problem related to time in quantum mechanics is the so-called "tunneling 
time problem". We can raise another, more general, question about the time. Let us 
consider a system which evolves in time. Let x is one of the observables of the system. 
During the evolution the value of x changes. We are considering a subset F of possible 
values of x- The question is how long the values of x belong to this subset. 
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There is another version of the question. If we know the final state of the system, we 
may ask how long the values of x belong to the subset under consideration when the system 
evolves from the initial to a final predefined state. The question about the tunneling time 
belongs to such class of the problems. Really, in the tunneling time problem we ask about 
the duration the particle spends in a specified region of the space, and we know that the 
particle has tunneled out, i.e., it is on the other side of the barrier. We can expect that 
such a question can not always be answered. Here our goal is to obtain the conditions 
under which it is possible to answer such question. 

One of the possibilities to solve the time problem is to answer what exactly the word 
"time" means. The meaning of every physical quantity is determined by the procedure of 
its measurement. Therefore, we have to construct a scheme of an experiment (this can be 
a gedanken experiment) which measures the quantity with the properties corresponding to 
the classical time. 

3.1 The model of the time measurement 

We consider a system evolving with time. Let one of the quantities describing the system to 
be X- Operator x corresponds to this quantity. For simplicity we assume that the operator 
X has a continuous spectrum. The case with discrete spectrum will be considered later. 

The measuring device interacts with the system only if x is in some region near the 
point xd, the concrete value of which depends on the detector only. If we want to measure 
the time when the system is in a large region of one has to use many detectors. In the 
case of tunneling a similar model was introduced by Steinberg [35] and developed in our 
paper [23]. The strong limit of such a model for analysis of the measurement effect for the 
quantum jumps has been used in Ref. [25]. 

We shall use the weak measurement concept described in Sec. El The operator A will 
be represented by operator 

^(Xd) = |xd)(xd| = 5(x-Xd). (8) 

It is assumed that after time t the readings of the detectors are collected and averaged. 

Hamiltonian ^ with A given by (jH)) represents a constant force acting on the detec- 
tor D when the quantity x is very close to xd- This force induces the change of the 
detector's momentum. From the classical point of view, the change of the momentum is 
proportional to the time the particle spends in the region around xd, and the coefficient 
of proportionality equals the force acting on the detector. We assume that the change of 
the mean momentum of the detector is proportional to the time the constant force acts on 
the detector and that the time the particle spends in the detector's region coincides with 
the time the force acts on the detector. 

We can replace the S function by the narrow rectangle of height 1/L and of width L in 
the X space. From Eq. Q it follows that the force acting on the detector when the particle 
is in the region around xd is = —X/L. The time the particle spends until time moment 
t in the unit-length region is 

r{t) = -UiP^it)) - M (9) 
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where (pq) and {pq{t)) are the mean initial momentum and the momentum after time t, 
respectively. If one wants to find the period the system spends in the region of the finite 
width, one must sum expressions type © many times. 

When the operator x has a discrete spectrum, one may ask how long the quantity x 
has the value xd- To answer this question the detector must interact with the system only 
when X = Xd- If this is satisfied, the operator D{xb) takes the form 

^(Xd) = |xd)(xd|. (10) 

The force, acting on the detector in this case equals to —A. The duration the quantity 
X has the value xd is given by Eq. too. Note that now formulae do not depend on 
spectrum of the operator x. 



3.2 The dwell time 

To shorten of the notation, the operator 

F{xB,t) = [ D{xB,h)dt, (11) 
Jo 

is introduced, where 

DixB,t) = Ulit)DixB)Usit). (12) 

After measurement, from the density matrix of the detector, in the first-order approx- 
imation we find that the average change of the detector momentum in the time interval t 
is — A(F(xd,^))- From Eq. (0) we obtain the dwell time until time moment t, 

T{x,t) = (^Fix,t)) . (13) 
Then the time spent in the region F is 

r(F;t) = / Tix,t)dx= f dt' [ dxP{x,t'), (14) 
Jr Jo Jr 

where Pix^t') = {D{x,t)) is the probability for the system to have the value x time 
moment t'. 

When X is the coordinate of the particle Eq. ()14j] yields the well-known expression for 
the dwell time [3,23]. This time is the average over the entire ensemble of the systems, 
regardless of their final states. 

A relationship between dwell, transmission, and reflection times recently has been anal- 
ysed in paper [36] while in paper [37] a relation between the group delay and the dwell time 
for quantum tunneling is derived. It is shown that the group delay is equal to the dwell 
time plus self-interference delay which depends on the dispersion outside the barrier. The 
analysis shows that there is nothing superluminal in quantum tunneling and the Hartman 
effect for tunneling quantum particles can be explained by the saturation of the integrated 
probability density under the barrier. 
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3.3 The definition of time under condition that the system is in 
the given final state 

Further, the case when the final state of the system is known will be considered . We may 
ask how long the values of x belong to the subset under consideration, T, on condition that 
the system evolves to the definite final state f. More specifically, we might know that the 
final state of the system belongs to a certain subspace TCf of system's Hilbert space. The 
projection operator that projects the vectors from the Hilbert space of the system into the 
subspace H{ of the final states will be denoted Pf. As far as the considered model gives the 
correct result for the time averaged over the entire ensemble of the systems, we can try to 
take the average only over the subensemble of the systems with the given final states. At 
first, the momenta Pq of each measuring device after the interaction with the system are 
measured. Subsequently, we perform the final, postselection measurement on the systems 
of the ensemble. Then we collect the outcomes Pq only of those systems the final state of 
which turns out to belong to the subspace Tif. 

Using Eq. ((Zj) in Sec. |2l we obtain the duration, on condition that the final state of the 
system belongs to the subspace Hf [22], 



nix^i) 



mit)) 

1 

+ 



Pfit)Fix,t) + Fix,t)m) 



((g)(pq)-Re(gpq))( Pf(t),F(x,t) 



(15) 



Eq. (fTH|l consists of two terms, thus, we can introduce two expressions with the dimen- 
sion of time 



T. 



(1) 



T, 



(2) 



ix,t) 
ix,t) 



mit)) 
1 



m),F{x,t) 



(16) 
(17) 



Then, the time the system spends in the subset F on condition that the final state of the 
system belongs to the subspace TCf can be rewritten in the form 



r^'\x,t) + I ((g) (pq) - Re{qp^))Tr{x,t). 



-(2), 



ri{x,t) 

The quantities (x, t) and (x, t) are related to the real and imaginary parts of the 
complex time introduced by Sokolovski et al [38]. In our model the quantity Tf (x, t) is real. 



contrary to the complex-time approach. The components of time t^^^ and 
too. Therefore, this time can be interpreted as the duration of an event. 

If the commutator Pf(t), F{x-, t) in Eq. ()15j) is not zero then, even in the limit of very 
weak measurement, the measured value depends on the particular detector used. This 
means that in such a case we cannot obtain a definite value for the conditional time. 
Moreover, the coefficient {{q){pq) — Re(gpq)) may be zero for the specific initial state of 
the detector, e.g., for the Gaussian distribution of the coordinate q and momentum p^. 



(2) 



are real, 
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The conditions to determine the time uniquely in a case when the final state of the 
system is known takes, thus, the form 



Pf(t),F(x,t)J =0 (19) 

which can be understood from on general principles of the quantum mechanics, too. Now, 
we ask how long the values of x belong to a certain subset when the system evolves to 
the given final state under assumption that the final state of the system is known with 
certainty. In addition, we want to have some information about the values of the quantity 
X- However, if the final state is known with certainty, we may not know the values of x iii 
the past and, vice versa, if we know something about we may not definitely determine 
the final state. Therefore, in such a case the question about the time when the system 
evolves to the given final state cannot be answered definitely and the conditional time has 
no reasonable meaning. 

The quantity rf(t) according to Eqs. (fT5|l and (fT6|l has many properties of the classical 
time. So, if the final states {/} constitute the full set, then the corresponding projection 
operators obey the equality of completeness J2fPf — ^- Then, from Eq. (|TH|l we obtain 
the expression 

j2{m)Mx,t) = Tix,t). (20) 
/ 

The quantity (-Pf(t)) is the probability that the system at the time t is in the state /. 
Eq. shows that the full duration equals the average over all possible final states, as it 
is a case in the classical physics. From Eq. (j20|l and Eqs. (fT6|) . f|T7jl it follows 

J2(m)rl'\x,t) = r{x,t), (21) 
/ 

J2im)^Pix,t) = 0. (22) 

/ 

We suppose that the quantities t^^^ {x, t) and r/^"* (x, t) can be useful even in the case when 
the time has no definite value, since in the tunneling time problem the quantities (fT?)|) and 
(I17|l correspond to real and imaginary parts of the complex time, respectively [23]. 

The eigenfunctions of the operator x constitute the full set / = 1; where the 

integral must be replaced by the sum for the discrete spectrum of the operator x- From 
Eqs. (jHl), (fTT|) . (fT3j) we obtain the equality 

rf(x,t)dx = t, (23) 

which shows that the time during which the quantity x has any value equals to t, as it is 
in the classical physics. 

3.4 Example: two-level system 

The obtained formalism can be applied to the tunneling time problem [23]. In this section, 
however, we will consider a simpler system than the tunneling particle, i.e., a two-level 
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system. The system is forced by the perturbation V that causes the jumps from one state 
to another. The time the system is in a given state will be calculated. 
The Hamiltonian of this system is 

H = Ho + V (24) 

where = ^(3"3/2 is the Hamiltonian of the unperturbed system and V = fO"+ + v*a- is 
the perturbation. Here ai, (T2, (T3 are Pauli matrices and a± = |((3"i±Z(T2). The Hamiltonian 
Hq has two eigenfunctions |0) and |1) with the eigenvalues —hjj/2 and hjj/2, respectively. 
The initial state of the system assumed to be |0). 

lYiom. Eq. ()13|1 we obtain the times the system spends in the energy levels and 1, 
respectively, 

^(1.*) = \i^-'€]t-^.Anmil-i^) (26) 



where f2 = uj"^ ^ 4(|t'|//i)2. From Eqs. (fTI)|) and (fTTj) we can obtain the conditional time. 
The components ()16|) and (jlTj) of the time the system spends in the level under condition 
that the final state after measurement is |1) are 

rf^(0,t) = \ (27) 

(2) - - ( ^ I VL 



When Vtt = 21111, where n G Z, the quantity t{ (0,t) tends to infinity. This happens 
because at these time moments the system is in the state |1) with the probability 0, and 
one cannot consider the interaction with the detector as very weak. 

On the other hand, the components of the time the system spends in level under 
condition that the final state is |0) are 

_ (1 + 3^) t + (1 - ^) sinm + t cosm) 
'° " 2((l + S) + (l-S)cos(^^t)) ' ^''^ 



(2), , . v^-n^)«in(f^)(^cos(ft)-^sin(ft)) 

^0 0,t = ^ -/ . ^ . 30 

2((l + S) + (l-f)cos(l]t)) 

The time the system spends in level 1 under condition that the final state is |0) may 
be expressed as 

+ t cos(r]t) - ^ sin(fit)) 

= 2((i + ^) + (i-^)cos(r^t)) • ^''^ 

The quantities T(0,t), T{l,t), T[^\o,t), TQ^\o,t) and TQ^\l,t) are shown in Fig.[T] The 
quantity Tq (0,t) is shown in Fig. |21 Note that the partial durations at the given final 
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1 2 3 4 5 6 

t 



Figure 1: The times the system spends in the energy level 0, T(0,t) (dashed line), and 
level 1, T(l,t) (dotted line), according to Eqs. and (|^. respectively. The quantity 
T[^\o,t), Eq. (P7|) . is shown as sohd straight line. The quantities TQ^\o,t) and TQ^\l,t), 
shown by curves 1 and 2, were calculated according to Eqs. and (jSH), respectively. 
The parameters are u = 2, Q = 4. 
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state are not necessarily monotonic as it is with the full duration, because the final state at 
different time moments can be reached by different paths. We can interpret the quantity 
TQ^\o,t) as the time the system spends in the level on condition that the final state is 
|0), but at certain time moments this quantity is greater than t. In such cases the quantity 
TQ^\l,t) becomes negative at certain times. This is a consequence of the fact that for 
the system under consideration the condition is not fulfilled. The peculiarities of the 
behavior of the conditional times show that it is impossible to decompose the unconditional 
time into two components having all classical properties of the time. 

4 Tunneling time 

The best-known problem of time in quantum mechanics is the so-called "tunneling time 
problem". This problem is still the subject of much controversy, since numerous theories 
contradict each other in their predictions for "the tunneling time" . Many of the theoretical 
approaches can be divided into three categories. First, one can study the evolution of the 
wave packets through the barrier and get the phase time. However, the correctness of 
the definition of this time is highly questionable [39]. Another approach is based on the 
determination of a set of dynamic paths, i.e., the calculation of the time the different 
paths spend in the barrier and averaging over the set of the paths. The paths can be 
found from the Feynman path integral formalism [38] , from the Bohm approach [40-45] , or 
from the Wigner distribution [46]. The third class uses a physical clock which is used for 
determination of the time elapsed during the tunneling (Biittiker and Landauer used an 
oscillatory barrier [39], Baz' suggested the Larmor time [47]). One more approach is based 
on a model for tunneling based on stochastic interpretation of quantum mechanics [48-51]. 

The problems rise also from the fact that the arrival time of a particle to a definite 
spatial point is a classical concept. Its quantum counterpart is problematic even for the 
free particle case. In classical mechanics, for the determination of the time the particle 
spends moving along a certain trajectory, one has to measure the position of the particle 
at two different moments of time. In quantum mechanics this procedure does not work. 
From Heisenberg's uncertainty principle it follows that we cannot measure the position of 
a particle without alteration of its momentum. To determine exactly the arrival time of 
a particle, one has to measure the position of the particle with great precision. Because 
of the measurement, the momentum of the particle will have a big uncertainty and the 
second measurement will be indefinite. If we want to ask about the time in quantum 
mechanics, we need to define the procedure of measurement. We can measure the position 
of the particle only with a finite precision and get a distribution of the possible positions. 
Applying such a measurement, we can expect to obtain not a single value of the traversal 
time but a distribution of times. 

In paper [52] the tunneling time distribution for photon tunneling is analysed theoret- 
ically as a space-time correlation phenomenon between the emission and absorption of a 
photon on the two sides of a barrier. The analysis is based on an appropriate counting 
rate formula derived at first order in the photon-detector interaction and used in treating 
space-time correlations between photons. 

There are two different but related questions connected with the tunneling time problem 
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[53]: 

(i) How much time does the tunnehng particle spend under the barrier? 

(ii) At what time does the particle arrive at the point behind the barrier? 

There have been many attempts to answer these questions. However, there are several 
papers showing that according to quantum mechanics the question (i) makes no sense 
[53-56]. Our goal is to investigate the possibility to determine the tunneling time using 
weak measurements. 



4.1 Determination of the tunneling time 



To answer the question of how much time does the tunneling particle spends under the 
barrier, we need a criterion of the tunneling. The following criterion is accepted: the 
particle had tunneled in the case when it was in front of the barrier at first and later it was 
found behind the barrier. We shall require that the mean energy of the particle and the 
energy uncertainty should be less than the height of the barrier. Following this criterion, 
the operator corresponding to the "tunneling-flag" observable is introduced 



MX)=Q{x-X), 



(32) 



where G is the Heaviside unit step function and X is a point behind the barrier. This 
operator projects the wave function onto the subspace of functions localized behind the 
barrier. The operator has two eigenvalues: and 1. The eigenvalue corresponds to 
the fact that the particle has not tunneled out, while the eigenvalue 1 corresponds to the 
appearance of particle behind the barrier. 

We will work with the Heisenberg representation. In this representation, the tunneling 
flag operator becomes 



fTit,X) = exp ( '-Htj /T(X)exp l-^Ht 



(33) 



To take into account all the tunneled particles, the limit t +oo must be taken. So, 
the "tunneling-flag" observable in the Heisenberg picture is represented by the operator 
/r(c)0,X) = limj^+oo fritjX). One can obtain the exphcit expression for this operator. 
The operator fT{t,X) obeys the standard equation 



zh^Jrit^X) 



fT{t,X),H 

The commutator in Eq. may be expressed as 



(34) 



fT{t,X),H 



exp I —Ht 

n 



fTiX),H 



exp 



h 



Ht 



If the Hamiltonian has the form H 



i^p^ + V{x), then the commutator becomes 



fT{X),H =ihJ{X), 



(35) 
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where j{X) is the probabihty flux operator, 



J{x) 



2M 



{x\p + p\x) {x\) . 



(36) 



Therefore, the following equation for the commutator can be written 

fT{t,X),Hj =thJ{X,t). (37) 

The initial condition for the function f{t,X) may be defined as 

/T(t = 0,X) = /T(X). 

From Eqs. ()34|1 and (j37jl we obtain the equation for the evolution of the tunneling-flag 
operator 

ih^jT{t,x)=thJ{x,t). 



(38) 

From Eq. ()38|1 and the initial condition, an explicit expression for the tunneling-flag oper- 
ator follows 



/r(t,X) = /T(X) 



J(X, ti)dti. 



(39) 



In the already mentioned question of how much time does the tunneling particle spend 
under the barrier, we shall be interested in those particles, which we know with certainty 
have tunneled out. In addition, we want to have some information about the location of 
the particle. However, one may ask whether the quantum mechanics allows one to have the 
information about the tunneling and location simultaneously? The projection operator 



DiT) 



\x) {x\dx 



(40) 



represents the probability for the particle to be in the region F. Here \x) is the eigenfunction 
of the coordinate operator. In the Heisenberg representation this operator takes the form 



^(F,^) = exp ( ^Htj Z)(F)exp (-^Ht 



(41) 



From Eqs. (^Hjl . (jSHI), and we see that the operators D(T,t) and /t(c)0,X) in general 
do not commute. This means that we cannot simultaneously have the information about 
the tunneling and location of the particle. If we know with certainty that the particle 
has tunneled out then we can say nothing about its location in the past, and if we know 
something about the location of the particle, we cannot determine definitely whether the 
particle has tunnel out. Therefore, the question of how much time does the tunneling 
particle spend under the barrier cannot have definite answer, if the question is so posed 
that its precise definition requires the existence of the joint probability that the particle 
is found in F at time t and whether or not it is found on the right side of the barrier at 
a sufficiently later time. A similar analysis has been performed in Ref. [56]. It has been 
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Figure 3: The configuration of tlie measurement of the tunnehng time. The particle P 
is tunnehng along x coordinate and weakly interacts with detectors D. The barrier is 
represented by the hatched rectangle. The interaction with the individual detectors occurs 
only in the narrow region limited by the horizontal lines. The changes of the momenta of 
the detectors are represented by arrows. 



shown that, due to noncommutability of operators, there exist no unique decomposition of 
the dwell time. 

This conclusion is, however, not negative altogether. We know that \x){x\dx = 1 

and 1, /t(oo, X) = 0. Therefore, if the region F is large enough, one has a possibility to 
answer the question about the tunneling time. 

/^From the fact that the operators -D(F, t) and /t(oo, X) do not commute we can predict 
that the measurement of the tunneling time will yield a value dependent on the particular 
detection scheme. We shall assume the detector is made so that it yields some value. 
But if we try to measure noncommuting observables, the measured values depend on the 
interaction between the detector and the measured system. So, in the definition of the 
Larmor time there is a dependence on the type of boundary attributed to the magnetic- 
field region [3]. 



4.2 The model of the time measurement 

We consider a model for the tunneling time measurement which is somewhat similar to the 
gedanken experiment used to obtain the Larmor time but is simpler and more transparent. 
This model was proposed by Steinberg [35] , however, it was treated in a nonstandard way, 
introducing complex probabilities. Here we shall use only the formalism of the standard 
quantum mechanics. 

Our system consists of a particle P and a number of detectors D [23]. Each detector 
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interacts with the particle only in the narrow region of space. The configuration of the 
system is shown in Fig. El When the interaction of the particle with the detectors is weak, 
the detectors do not infiuence the state of the particle. Therefore, we can analyze the 
action of the detectors separately. This model is a particular case of time measurement 
presented in Sec. with xd being the position of the detector xd- Similar calculations 
were done for detector's position rather than momentum by lannaccone [57]. 

In the moment t = the particle is before the barrier, therefore, {x\pp{0)\x') 7^ only 
when a; < and x' < 0, where pp(0) is the density matrix of the particle P. 



4.3 Measurement of the dwell time 

As in Sec. 13.21 we obtain the time the particle spends in the unit length region between 
time instances and t 

r°-(x,t) = (F(a;,t)). (42) 
The time spent in the space region restricted by the coordinates xi and X2 is 

PX^ PX2 PCO 

t°™(x2,xi)= / r°^(x,t ^ oo)dx = dx p{x,t)dt (43) 

J xi J xi Jo 

which is a well-known expression for the dwell time [3]: the dwell time is the average over 
an entire ensemble of particles regardless whether they tunneled or not. The expression for 
the dwell time obtained in our model is the same as the well-known expression obtained by 
other authors. Therefore, we can expect that our model can yield a reasonable expression 
for the tunneling time as well. 



4.4 Conditional probabilities and the tunneling time 

Having seen that our model is capable to give the time averaged over entire ensemble of 
the particles, let us now take the average over the subensemble of the tunneled particles 
only. This will be done similarly to Sec. l3.3l with Pf replaced by the tunneling-fiag operator 
/t(X) defined by Eq. From Eq. (fT3|) we obtain the duration the tunneled particle 

spends in the unit length region around x until time t [23] 



= :;7T7r^{Mt,X)F{x,t) + F{x,t)Mt,X) 

+Tr7T^-— 7 m iP,) - Mqp,)) ( iMt, X), Fix, t)] ) . (44) 

The obtained expression for the tunneling time is real, contrary to the complex-time 
approach. It should be noted that this expression, even in the limit of weak measurement, 
depends on a particular detector. If the commutator [fT(t, X), F{x,t)] is zero, the time 
has a well-defined value. If the commutator is not zero, only the integral of this expression 
over a large region has meaning of an asymptotic time related to the large region as we 
will see in Sec. 14.71 



15 



Equation can be rewritten as a sum of two terms, the first term being independent 
and the second dependent on the detector, i.e.. 



T{x,t) = rT-(x,t) + - ((g) (p,) - Re{qp,))T^^{x,t) (45) 

where 

= o/.- u (fTit,X)Fix,t)+Fix,t)Mt,X)), (46) 

r^^{x,t) = (\fT{t,X),F{x,t)]) . (47) 

The quantities r'^™(x, t) and T]-^f{x,t) are independent of the detector. 

In order to separate the tunneled and reflected particles the limit t ^ oo should be 
taken. Otherwise, the particles that tunneled after the time t will not contribute. If we 
introduce the operators 



F{x) = / D{x,h)dh, (48) 

^0 

POO 

N{x) = / J{x,ti)dti. (49) 



then from Eq. follows that the operator /j'(oo,X) is friX) + N{X). If the particle 
before the barrier is initially, then 

/t(X)/)p(0) = /)p(0)/t(X) = 0. 

In the limit t — oo tunneling times become 

r^-(x) = — ^(iV(X)F(x)+F(x)iV(X)), (50) 



2{N{X)) 
1 / 

2t{NixY) 



^^'^l^) = :rT7^777Z7{N{X),F{x) ). (51) 



Let us define an "asymptotic time" as the integral of r(x, oo) over a wide region con- 
taining the barrier. Since the integral of T^rr (a;) is very small compared to that of r"'""'^(x) 
as we shall see later, the asymptotic time is effectively the integral of t'^'^^Ix) only. This 
allows us to identify r^"'^(x) as "the density of the tunneling time". 

In many cases for the simplification of mathematics it is common to write the integrals 
over time as the integrals from — oo to +oo. In our model we cannot, without additional 
assumptions, integrate Eqs. (pSj) . (jl^ from — oo because the negative time means the 
motion of the particle to the initial position. If some particle in the initial wave packet had 
negative momentum then in the limit t — >^ — oo it was behind the barrier and contributed 
to the tunneling time. 
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4.5 Properties of the tunneling time 

As stated, the question of how much time does a tunneling particle spend under the barrier 
has no exact answer. We can determine only the time the tunneling particle spends in a 
large region containing the barrier. In our model this time is expressed as an integral of 
quantity (jMIj) over this region. In order to determine the properties of this integral it is 
useful to determine the properties of the integrand. 

To be able to expand the range of integration over time to — oo, it is necessary to have 
the initial wave packet far to the left from the points under the investigation and this wave 
packet must consist only of the waves moving in the positive direction. 

It is convenient to perform calculations in the energy representation. Eigenf unctions 
of the Hamiltonian Hp are \E,a), where a = ±1. The sign '+' or '— ' corresponds to the 
positive or negative initial direction of the wave, respectively. Outside the barrier these 
eigenfunctions are 



M 



x\E,+) = { ; - - ; " (52) 



(exp (jiPex) + r{E) exp [-\pex)) , x < 0, 



M 



M 



t{E) exp [j:pEx) , x>L 



. ,2.npJiE)exp{-lpEx), x<0, 

^ (exp i-lpEx) - ^f\E) exp , x>L 

where t{E) and r{E) are transmission and reflection amplitudes respectively, and 



Pe = V2ME. (54) 

M is the mass of the particle. The barrier is in the region between x = and x = L. 
These eigenfunctions are orthonormal, i.e., 

{E,a\E',a') = 6a,a'S{E - E'). (55) 

The evolution operator is 

Up{t) = ^j^ \E,a){E,a\eyip(^^-E?jdE. 

Then the operator F{x) assumes the form 

F{x) = j dti^jjdE dE'\E, a) {E, a\x) {x\E', a') {E', a'\ exp (^^ {E - E') t^ 

where the integral over the time yields 2Trh6{E — E') and, therefore, 
F{x) = 27rn^ j dE\E,a){E,a\x){x\E,a'){E,a'\. 



Similarly, we find 



N{x) = 27thY^ j dE\E,a){E,a\J{x)\E,a'){E,a'\ 
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If the initial wave packet consisting only of the waves moving in the positive direction is 
assumed, then one has 



(iV(x)) = 27ihjdE(^\E,+){E,+\J{x)\E,+){E,+\ 



From the condition X > L it follows that 



(iV(X)) = I dE{\E,+)\m\'{E,+\)- (56) 
For X < we obtain the following expressions for the quantities T^'^^^Xjt) and Tj^'^{x,t) 



+ r*{E)exp y^^PExjj {E,+\j, (57) 
r^M ^ ^/dE(|E,+)^|t(E)p(.(E)exp(-2lp.. 



- r*{E)exp[2-pEx] ] {E,+\). (5J 



For X > L these expressions take the form 

T'^^^(x,t) = I AeI — |t(^)P (^2- ^^r*(E)exp ( 



^ ^^^r(E)exp (-2'-pex ) ) ), (59) 



t(E) " V ^ 

rT-(a;,t) = / dEl\E,+) — \t{E)\' ( ^r*{E)exp ( 2^p^x') 

2{N{X))J \ Pe \t*iE) "^yh"^" J 

^*^^KiE)exp (-2lpEx]) {E,+\). (60) 



To illustrate the obtained formulae, the 5-function barrier 

V{x) = n6{x) 

and the rectangular barrier will be used. The Gaussian incident wave packet initially is far 
to the left of the barrier. 

In Fig. m and we see interferencelike oscillations near the barrier. Oscillations are 
present not only in the front of the barrier but also behind the barrier. When x is far 
from the barrier the "time density" tends to a value close to 1. This is in agreement 
with classical mechanics because in the chosen units the mean velocity of the particle is 
1. Fig. El shows additional property of "tunneling time density": it is almost zero in the 
barrier region. This explains the Hartmann and Fletcher effect [58,59]: for opaque barriers 
the effective tunneling velocity is very large. 



18 



2.0 




0.0 I I 

-15 -10 -5 5 10 15 

X 



Figure 4: The asymptotic time density for (5-function barrier with the parameter Q = 2. 
The barrier is located at the point x = 0. The units are such that h = 1 and M = 1 and the 
average momentum of the Gaussian wave packet (p) = 1. In these units length and time 
are dimensionless. The width of the wave packet in the momentum space is u = 0.001. 
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Figure 5: The asymptotic time density for rectangular barrier. The barrier is localized 
between the points x = and x = 5 and the height of the barrier is Vq = 2. The used 
units and parameters of the initial wave packet are the same as in Fig. HI 
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4.6 The reflection time 

We can easily adapt our model for the reflection too. In doing this, one should replace the 
tunneling-flag operator Jt by the reflection-flag operator 

fR = l- It. (61) 
Replacement of /t by fR in Eqs. (j50|) and (j5T| gives 

Cfa{t = oo,X))r^^fl(x) = r^^^ix) - {Mt = oo, X))r^-\x). (62) 
We see that in our model the important condition 

where T and R are the transmission and reflection probabilities is satisfied automatically. 

If the wave packet consists of waves moving in the positive direction, the density of 
dwell time becomes 

r^'"{x,t) = 27ih j dE{\E,+){E,+\x){x\E,+){E,+\). (64) 

For X < we have 

rD"(a;,t) = M j dE (\E,+)^(^ + \r{E)\'^ + r{E)e^^(^2^PEX^ 

+ r*(E)exp (^'-pE^^ (65) 

and for the reflection time we obtain the "time density" 

M 

1 - {N{X)) 

+ ^ (1 + |r(E)|2) r{E) exp (^-2^psx^ 



f dEl\E,+)—{2\r{E)f 
J \ Pe 



+ r*{E)exp { 2^peX ] ] {E,+\) . (66) 



For X > L the density of the dwell time is 

r°"(a;,t) = M [ dE ( \E, +) — \t{E)\^{E, +\ 
J \ Pe 

and the "density of the reflection time" may be expressed as 



(67) 



n^) = f |dE^|E,+)^|t(i?)p(^r*(i?)exp(2lp^x 



+ yir^(^) ( -2^Pi?x ) ) (E, +1 ) . (68) 
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Figure 8: Reflection time density for rectangular barrier in the region behind the barrier. 
The parameters and the initial conditions are the same as in Fig. 

We will illustrate the properties of the reflection time for the same barriers and Gaussian 
incident wave packet initially localized far to the left from the barrier. In Figs. IHl and 
[7J one can see the interference-like oscillations at both sides of the barrier. Since for 
the rectangular barrier the "time density" behind the barrier is very small, this part is 
presented in Fig. |H1 Behind the barrier, the "time density" at certain points becomes 
negative. This is because the quantity t^°^{x) is not positive definite. Nonpositivity is the 
direct consequence of noncommutativity of the operators in Eqs. ()50|) and ()51|) . There is 
nothing strange in the negativity of r^'^^(a;) because this quantity has no physical meaning. 
Only the integral over the large region has the meaning of time. When x is far to the left 
from the barrier the "time density" tends to a value close to 2 and when x is far to the 
right from the barrier the "time density" tends to 0. This is in agreement with classical 
mechanics because in the chosen units, the velocity of the particle is 1 and the reflected 
particle crosses the area before the barrier two times. 

4.7 The asymptotic time 

As mentioned above, we can determine only the time that the tunneling particle spends 
in a large region containing the barrier, i.e., the asymptotic time. In our model this time 
is expressed as an integral of quantity (j50|) over this region. We can do the integration 
explicitly. 

The continuity equation yields 



d 



J{xD,t) = 0. 



(69) 



dxD 



The integration can be performed by parts 
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If the density matrix pp(0) represents localized particle then linit^oo t)pp(O) j = 0. 

Therefore we can write an effective equality 

/•OO _ Q l>00 

J D{xD,ti)dh = — j tiJ{xD,h)dh. (70) 
We introduce the operator 

/■oo 

f{x) = / tiJ{x,ti)dti. (71) 



We consider the asymptotic time, i.e., the time the particle spends between points xi and 
X2 when xi —oo, X2 — > +oo, 

t^"°(a;2,a;i) = f \^''''{x)dx. 

J Xl 

After the integration we have 

t'^"(x2, Xl) = tT"'^(x2) - t^'^'^(Xi) (72) 

where 

t^^'^'i^) = ^^^^ (iV(x)f (x) + f (x)iV(x)) . (73) 

If we assume that the initial wave packet is far to the left from the points under the 
investigation and consists only of the waves moving in the positive direction, then Eq. ()72|1 
may be simplified. 

In the energy representation the operator (f7T|) is 

f{x)= J hdhJ^JJ dEdE'\E,a){E,a\J{x)\E',a){E',a'\exp(^^{E - E')ti^ . 



The integration over time yields 2i'Kh'^-^5{E — E') and we obtain 
f{x) = -ih27ThJ2 J '^E\E,a) (-^{E,a\J{x)\E',a') 

a, a' 



(iV(X)f(x)) 



{E,a'\ 

E'=E 



+ {E,a\J{x)\E,a')-^{E,a'\ 
-iM7r^h'^J2 j dE{^>\E,+){E,+\J{X)\E,a) 



X 



a 

^ -{E,a\J{x)\E',+) 



dE 



^^_^ + {E,a\J{x)\E,+)^] {E,+\m). 



Substituting expressions for the matrix elements of the probability flux operator we obtain 
equation 

(iV(X)f(x)) = j dEmE,+)e{E)-^^t{E){E,+\^>) 

+Mx I dE{m\E,+) — \t{E)\^{E,+\^!) 
J Pe 

+ih^ j dE{^E,+)^r*{E)t\E)exp (^^^PeX^ (^,+1^)- 
23 



M 



When X ^ +00, the last term vanishes and we have 

{N{X)f{x)) = j dEmE,+)f{E)-^^t{E){E,+\^>) 

+Mx ! AEi^<\E,+) — \t{E)\^iE,+m, x^+00. (74) 
J Pe 

This expression is equal to (T(x)), 

(iV(X)f (x)) ^ {f{x)), X +00. (75) 
When the point with coordinate x is in front of the barrier, expression ()74|1 becomes 

(iV(X)f(x)) = -th [dE{<il\E,+)\t{E)\''(^—x 

J \npE 

iE) exp (-;^2psx^ + {E, +|^). 
When |x| is large the second term vanishes and we have 

(iV(X)f(x)) ^ Mx [ dE{^\E,+) — \t{E)mE,+\^) 

J Pe 

+ 1 dE{<il\E,+)\tiE)\'^^{E,+\<i^). (76) 

The imaginary part of expression (f7B|) is not zero. This means that for determination of 
the asymptotic time it is insufficient to integrate only in the region containing the barrier. 
For quasimonochromatic wave packets from Eqs. (|7T|l . (|7^ . ((ZBI), and (ffH|l we obtain 
the limits 

t^""(x2,xi) ^ 4'' + —M{x2-xi), (77) 

Pe 

C"(^2,a:i) ^ -t^- (78) 

where 



is the phase time and 



t^T- = h^{aTgt{E)) (79) 

t'^ = h^{\n\t{E)\) (80) 

is the imaginary part of the complex time. 

In order to take the limit a; —00 we have to perform more accurate calculations. 
The range of integration over time to cannot be extended —00 because such extension 
corresponds to the initial wave packet being infinitely far from the barrier. We can extend 
the range of the integration over the time to —00 only in N{X). For a; < we obtain the 
following equation 

(iv(x)f(x)) = ^J^m (^i*(x,t)^hix,t)-hix,t)^r,ix,t)^ (81) 

24 



where 



h{x,t) = J dE^\tiE)\^exp(^^{pEX- Et)^ {E,+\<i/), (82) 
l2ix,t) = j ^E^= (^exp {jPeX^ +r(E)exp (^^-pEX^^ exp [-j^E^j {E,+\%p) 

Ii{x,t) is equal to the wave function at the point x and the time moment t, when the 
propagation is in the free space and the initial wave function in the energy representation 
is |t(E)p(E, +1^). When t > and x -oo, then Ii{x,t) 0. That is why the 
initial wave packet contains only the waves moving in the positive direction. Therefore 
{N{X)T{x)) when x — > — oo. From this analysis it follows that the region in which 
the asymptotic time is well determined has to include not only the barrier but also the 
initial wave packet region. 

In such a case from Eqs. (fT2|) and (|T3| we obtain expression for the asymptotic time 



dE{<f\E, +)t*{E) (^^X2 - t{E){E, +1^) 



t'^"(a;2,xi^-oo) = -^— — I dE{<f\E,+)t*{E)[—X2-ih—]tiE){E,+\-^). (84) 

(TV (A)) 

From Eq. ()75p it follows that 

tT-(a;2,a:i -oo) = --^(^(x^)) (85) 

(7V(A j) 

where T(x2) is defined as the probability flux integral (|7T|l . Equations (jHH) and (fHHjl give 
the same value for tunneling time as does an approach in Refs. [60,61] 

The integral of quantity Tj^f{x) over a large region is zero. We have seen that it is not 
enough to choose the region around the barrier — this region has to include also the initial 
wave packet location. This fact will be illustrated by numerical calculations. 

The quantity T^f{x) for (5-function barrier is represented in Fig. El We see that t^^^{x) 
is not equal to zero not only in the region around the barrier but also it is not zero in the 
location of the initial wave packet. For comparison, the quantity r'^'^°(a;) for the same 
conditions is represented in Fig. ^1 



5 Arrival time 

The detection of the particles in time-of-fiight and coincidence experiments are common, 
and quantum mechanics should give a method for the calculation of the arrival time. The 
arrival time distribution may be useful in solving the tunneling time problem, as well. 
Therefore, the quantum description of arrival time has attracted much attention [62-77]. 
Aharonov and Bohm introduced the arrival time operator [62] 

Tab = ^(iX- x)- + -{X - x)] . (86) 
2 V P P J 

By imposing several conditions (normalization, positivity, minimum variance, and sym- 
metry with respect to the arrival point X) a quantum arrival time distribution for a free 



25 




Figure 9: The quantity tJ™^{x) for 6 function barrier with the parameters and initial 
conditions as in Fig. HI The initial packet is shown by dashed line. 




Figure 10: Tunneling time density for the same conditions and parameters as in Fig. 121 
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particle was obtained by Kijowski [63]. Kijowski's distribution may be associated with 
the positive operator valued measure generated by the eigenstates of Tab- However, Ki- 
jowski's set of conditions cannot be applied in a general case [63]. Nevertheless, arrival 
time operators can be constructed even if the particle is not free [73,78]. 

Since the mean arrival time even in classical mechanics can be infinite or the particle 
may not arrive at all, it is convenient to deal not with the mean arrival time and corre- 
sponding operator T, but with the probability distribution of the arrival times [24]. The 
probability distribution of the arrival times can be obtained from a suitable classical defi- 
nition. The noncommutativity of the operators in quantum mechanics is circumvented by 
using the concept of weak measurements. 

5.1 Arrival time in classical mechanics 

In classical mechanics the particle moves along the trajectory H{x, p) = const as t increases. 
This allows us to work out the time of arrival at the point x{t) = X, by identifying the 
point (xo,Po) of the phase space where the particle is at t = 0, and then following the 
trajectory that passes by this point, up to arrival at the point X. If multiple crossings are 
possible, one may define a distribution of arrival times with contributions from all crossings, 
when no distinction is made between first, second and nth arrivals. In this article we will 
consider such a distribution. 

We can ask whether there is a definition of the arrival time that is valid in both classical 
and quantum mechanics. In our opinion, the words "the particle arrives from the left at 
the point X at the time t" mean that: (i) at time t the particle was in the region x < X 
and (ii) at time t + At {At 0) the particle is found in the region x > X. Now we apply 
the definition given by (i) and (ii) to the time of arrival in the classical case. 

Since quantum mechanics deals with probabilities, it is convenient to use probabilistic 
description of the classical mechanics, as well. Therefore, we will consider an ensemble of 
noninteracting classical particles. The probability density in the phase space is p{x,p;t). 

Let us denote the region a; < X as Fi and the region x > X as The probability 
that the particle arrives from region Fi to region r2 at a time between t and t + At is 
proportional to the probability that the particle is in region Fi at time t and in region F2 
at time t + At. This probability is 

n+(t)At = -^ / dpdxp{x,p;t), (87) 

where is the constant of normalization and the region of phase space Q has the following 
properties: (i) the coordinates of the points in Q are in the space region Fi and (ii) if the 
phase trajectory goes through a point of the region Q at time t then the particle at time 
t + At is in the space region F2. Since At is infinitesimal, the change of coordinate during 
the time interval At is equal to ;^At. Therefore, the particle arrives from region Fi to 
region F2 only if the momentum of the particle at the point X is positive. The phase space 
region Q consists of the points with positive momentum p and with coordinates between 
X — p/mAt and X. Then from Eq. (jHTj) we have the probability of arrival time 

n+(t)At = — / dp dxp{x,p-t). (88) 
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Since At is infinitesimal and the momentum of every particle is finite, we can replace x in 
Eq. (jHH|) by X and obtain the equality 



n+(t,X) = -^/ ^p{X,p;t)dp. (89) 
Jo ^ 

The obtained arrival time distribution n_|_(t, X) is well known and has appeared quite often 
in the literature (see, e.g., the review [73] and references therein). 
The probability current in classical mechanics is 

J{x-t) = [ ^p{x,p;t)dp. (90) 

From Eqs. ()89j) and (jHUIl it is clear that the time of arrival is related to the probability 
current. This relation, however, is not straightforward. We can introduce the "positive 
probability current" 

P 

J+{x]t)= I —p{x,p;t)dp (91) 
Jo m 

and rewrite Eq. (IHUj) as 

n+(t,X) = ^J+(X;t). (92) 

The proposed [79-81] various quantum versions of J+ even in the case of the free particle 
can be negative (the so-called backflow effect). Therefore, the classical expression ()92|) for 
the time of arrival becomes problematic in quantum mechanics. 

Similarly, for arrival from the right we obtain the probability density 

n_(t,X) = ^J„(X;t), (93) 

where the negative probability current is 

J_(x;t)= / ^p(x,p;t)dp. (94) 

We see that our definition given at the beginning of this section leads to the proper 
result in classical mechanics. The conditions (i) and (ii) does not involve the concept of 
the trajectories. We can try to use this definition also in quantum mechanics. 



5.2 Weak measurement of arrival time 

The proposed definition of the arrival time probability distribution can be used in quantum 
mechanics only if the determination of the region in which the particle is does not disturb 
the motion of the particle. This can be achieved using the weak measurements of Aharonov, 
Albert and Vaidman [16-21]. 

We use the weak measurement, described in Sec. |2] The detector interacts with the 
particle only in region Fi. As regards the operator A we take the projection operator Pi 
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which projects into region Fi. In analogy to Ref. [17], we define the "weak value" of the 
probability of finding the particle in the region Fi , 



W{1) = (A) = ^^"^ (95) 

AT 

In order to obtain the arrival time probability using the definition from Sec. 15.11 we 
measure the momenta Pq of each detector after the interaction with the particle. After 
time At we perform the final, postselection measurement on the particles of our ensemble 
and measure if the particle is found in region F2. Then we collect the outcomes Pg only for 
the particles found in region F2. 

The projection operator projecting into the region F2 is P2. In the Heisenberg repre- 
sentation this operator is 

P2W = U{t)^p2Uit), (96) 

where U is the evolution operator of the free particle. Taking the operator B from Sec. El 
as P2{At) and using Eq. we can introduce a weak value H^(l|2) of probability to find 
the particle in the region Fi on condition that the particle after time At is in the region 
F2. The probability that the particle is in region Fi and after time At it is in region F2 
then equals 

W{1,2) = W{2)W{1\2). (97) 

When the measurement time r is sufficiently small, the inffuence of the Hamiltonian of the 
particle can be neglected. Using Eq. from Sec. El we obtain 

W{1,2) ^ ^(P2(At)A + AA(At)) + ^ {{pg){q) - Re{qpq)) ([A, A(At)]). (98) 

The probability 1^(1,2) is constructed using conditions (i) and (ii) from Sec. 15.11 the 
weak measurement is performed to determine if the particle is in the region Fi and after 
time At the strong measurement determines if the particle is in the region F2. Therefore, 
according to Sec. 15.11 the quantity W{1,2) after normalization can be considered as the 
weak value of the arrival time probability distribution. 

Equation ()98|1 consists of two terms and we accordingly can introduce two quantities 



and 



n^'^ = ^(AA(At) + A(At)A) (99) 



n('^ = ^([A,A(At)]). (100) 



Then 

2At 

W{1, 2) = n«At - -J- {{pg){q) - Re{qpq)) U^'l (101) 

If the commutator [A, A (At)] in Eqs. (jnn|) - p01|l is not zero, then, even in the limit 
of the very weak measurement, the measured value depends on the particular detector. 
This fact means that in such a case we cannot obtain a definite value for the arrival time 
probability. Moreover, the coefficient {{pq){q) — Re{qpq)) may be zero for a specific initial 
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state of the detector, e.g., for a Gaussian distribution of the coordinate q and momentum 

Pq- 

The quantities W(l,2), U^^^ and H^^) are real. However, it is convenient to consider 
the complex quantity 

Uc = +zn(2) = i^(AA(Ai)). (102) 

We call it the "complex arrival probability" . We can introduce the corresponding operator 

n+ = ^AP2(At). (103) 

By analogy, the operator 

ft- = ^AA(Ai) (104) 

corresponds to arrival from the right. 

The introduced operator n+ has some of the properties of the classical positive proba- 
bility current. Prom the conditions A + -A = 1 and A(At) + A(Ai) = 1 we have 

n+-n_ = i^(A(At)-A). 

In the limit — > we obtain the probability current J = hmAt^o(n+ — ft_), as in classical 
mechanics. However, the quantity (n+) is complex and the real part can be negative, in 
contrast to the classical quantity J+. The reason for this is the noncommutativity of the 
operators A and P2{At). When the imaginary part is small, the quantity (11+) after 
normalization can be considered as the approximate probability distribution of the arrival 
time. 

5.3 Arrival time probability 

The operator n_|_ was obtained without specification of the Hamiltonian of the particle and 
is suitable for free particles and for particles subjected to an external potential as well. In 
this section we consider the arrival time of the free particle. 

The calculation of the "weak arrival time distribution" W{1,2) involves the average 
value (11+). Therefore, it is useful to have the matrix elements of the operator 11+. It 
should be noted that the matrix elements of the operator 11+ , as well as the operator itself, 
are only auxiliary quantities and do not have an independent meaning. 

In the basis of momentum eigenstates \p), normalized according to the condition (pi|p2) = 
27rM(pi — P2), the matrix elements of the operator 11+ are 

(Pi|n+|p2) = -J-(pi|At>(At)tAt/(At)|p2) 



At 



1 pX poo 2 

= — / dxi / dx2e-sPi"i(a;i|t/(At)t|a;2)e^*'^^^-^^^*. (105) 
At J_oo Jx 



After performing the integration one obtains 

(Pi|n+|p2) = ^A^(^r^-p(^(^^^-^^0^ 
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Figure 11: The real part of the quantity {p\fl^\p), according to Eq. ()108|) . The correspond- 
ing classical positive probability current is shown with the dashed line. The parameters 
used are h = 1, m = 1, and At = 1. In this system of units, the momentum p is dimen- 
sionless. 



e*««-i.e*|-p,,/||i -e* (106) 



where y/i = exp(i7r/4). When 



the matrix elements of the operator 11+ are 

{p,\UM ^ ^^e^P . (107) 

This equation coincides with the expression for the matrix elements of the probability 
current operator. 

^From Eq. p06p we obtain the diagonal matrix elements of the operator n_|_ , 



(p\Il+\p) = — erfc + e"^^^*. (108) 

The real part of the quantity (p|n+|p) is shown in Fig.^^and the imaginary part in Fig.lT^ 

Using the asymptotic expressions for the error function erfc we obtain from Eq. 
that 

p 

lim (p|n+|p) — > — 
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Figure 12: The imaginary part of the quantity (p|n_|.|p). The parameters used are the 
same as in Fig. ITT] 

and (p|n+|p) 0, when p — oo, i.e., the imaginary part tends to zero and the real 
part approaches the corresponding classical value as the modulus of the momentum \p\ 
increases. Such behaviour is evident from Figs. II II and 1121 also. 

The asymptotic expressions for function erfc are valid when the argument of the erfc is 

large, i.e., \p\^f^ > 1 or 

At > A. (109) 

Here is the kinetic energy of the particle. The dependence of the quantity Re(p|n+|p) 
on At is shown in Fig. ^| For small At the quantity (p|n_|_|p) is proportional to l/v^At. 
Therefore, unlike in classical mechanics, in quantum mechanics At cannot be zero. Equa- 
tion imposes the lower bound on the resolution time At. It follows that our model 
does not permit determination of the arrival time with resolution greater than h/Ek- A 
relation similar to Eq. ()109p based on measurement models was obtained by Aharonov et 
al. [30] . The time-energy uncertainty relations associated with the time of arrival distribu- 
tion are also discussed in Refs. [69, 82] 



6 Summary 

The review and generalization of the theoretical analysis of the time problem in quantum 
mechanics and weak measurements are presented. The tunneling time problem is a part of 
this more general problem. The problem of time is solved adapting the weak measurement 
theory to the measurement of time. In this model the expression ()13|1 for the duration, when 
the arbitrary observable x has the certain value, is obtained. This result is in agreement 
with the known results for the dwell time in the tunneling time problem. 

Further we consider the problem of the duration when the observable x has a certain 



32 



1.5 




0.5 r I 

1 2 3 4 5 6 

M 

Figure 13: The dependence of the quantity Re(p|n+|p) according to Eq. ()108|) on the 
resolution time At. The corresponding classical positive probability current is shown with 
the dashed line. The parameters used are ^ = 1, m = 1, and p = 1. In these units, the 
time At is dimensionless. 

value on condition that the system is in the given final state. Our model of measurement 
allows us to obtain the expression (fTH|l of this duration as well. This expression has many 
properties of the corresponding classical time. However, such a duration not always has 
the reasonable meaning. It is possible to obtain the duration the quantity x has the certain 
value on condition that the system is in a given final state only when the condition ()19|) is 
fulfilled. In the opposite case, there is a dependence in the outcome of the measurements 
on particular detector even in an ideal case and, therefore, it is impossible to obtain the 
definite value of the duration. When the condition (fT^ is not fulfilled, we introduce two 
quantities (fT^ and (fTTj). characterizing the conditional time. These quantities are useful 
in the case of tunnehng and we suppose that they can be useful also for other problems. 

In order to investigate the tunneling time problem, we consider a procedure of time 
measurement, proposed by Steinberg [35]. This procedure shows clearly the consequences 
of noncommutativity of the operators and the possibility of determination of the asymptotic 
time. Our model also reveals the Hartmann and Fletcher effect, i.e., for opaque barriers 
the effective velocity is very large because the contribution of the barrier region to the time 
is almost zero. We cannot determine whether this velocity can be larger than c because 
for this purpose one has to use a relativistic equation (e.g., the Dirac equation). 

The definition of density of one sided arrivals is proposed. This definition is extended to 
quantum mechanics, using the concept of weak measurements by Aharonov et al [16-21]. 
The proposed procedure is suitable for free particles and for the particles subjected to an 
external potential, as well. It gives not only a mathematical expression for the arrival time 
probability distribution but also a way of measuring the quantity obtained. However, this 
procedure gives no unique expression for the arrival time probability distribution. 

In analogy with the complex tunneling time, the complex arrival time "probability 
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distribution" is introduced (Eq. ()102p ). It is shown that the proposed approach imposes 
an inherent hmitation, Eq. ()109|) . on the resolution time of the arrival time determination. 
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